Supersolid phase in spin dimer XXZ systems under magnetic field 
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Using quantum Monte Carlo method, we study, under external magnetic fields, the ground state 
phase diagram of the two-dimensional spin 5= 1/2 dimer model with an anisotropic intra-plane 
antiferromagnetic coupling. With the anisotropy 4 > A > 3, a supersolid phase characterized 
by a non-uniform bose condensate density that breaks translational symmetry is found. The rich 
phase diagram also contains a checkerboard solid and two different types of superfluid phase formed 
hj Sz = -1-1 and S'z = spin triplets, with finite staggered magnetization in z-axis and in-plane 
direction, respectively. As we show, the model can be realized as a consequence of including the 
next nearest neighbor coupling among dimers and our results suggest that spin dimer systems may 
be an ideal model system to study the supersolid phase. 

PACS numbers: 75.10.Jm, 75.45.+J, 05.30.Jp, 75.40.Mg 
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Five decades ago, Penrose and Onsager 'l| showed 
that the coexistence of both off-diagonal long-range or- 
der (ODLRO) and diagonal long range order (DLRO), i.e. 
a supersolid (SS) state, is impossible in commensurate 
solids. However, it was later demonstrated |2| that vacan- 
cies in solid may still exhibit Bose-Einstein condensation. 
And more recently, Prokof'ev and Svistunov 'i|| provide 
a proof that the presence of vacancies and/or interstitial 
atoms is a necessary condition for a SS. They therefore 
exclude the possibility of a bulk SS in the recent experi- 
mental observations of nonclassical moment of inertia in 
solid "^He |j|. Instead, interpretation based on superfluid 
interfaces is suggested. While experimental evidence of 
SS is still inconclusive, numerical stimulations have al- 
ready confirmed the presence of SS of bosons in various 
conditions, that is including next nearest neighbor inter- 
actions (n.n.n.) [5| or with softcore |g] in square lattices, 
and with hardcore in triangular lattices [3l- 

All of the models discussed above concentrate on 
studying bosons. Knowing the possibility of representing 
bosons by spins or vice versa, we wonder if SS phase also 
exists in spin systems. Recent experiments on spin dimer 
compounds, such as TlCuCla, KCuCle and BaSiCu206 
[81, have confirmed the formation of super fluid (SF) of 
the spin triplets under an external magnetic field. These 
spin bosons have also previously been found to be crys- 
tallized, signaled by the magnetization plateau, in dimer 
compounds such as SrCu2(B03)2^9']. The competition 
between repulsive interaction and kinetic motion decides 
whether the ground state is crystalized or is a SF llOJ. 
This also raises the possibility of a SS state. While SS is 
unstable against phase separation for hardcore boson in 
a square lattice la, lll| , the spin bosons are shown to be 
semi- hardcore |l2| . In this work, using quantum Monte 
Carlo (QMC) method we demonstrate that the ground 
state of a spin dimer XXZ model could be SF of two 
kinds of triplets, a quantum solid or, most importantly, 
a SS. This XXZ anisotropy can be easily realized if the 
n.n.n. interaction among dimers is considered. There- 



fore, while previous proposal of detecting SS in optical 
lattices [3 is technically difficult, we suggest that, due to 
the semi-hardcore nature of spin bosons, it may be more 
natural to look for a SS in spin dimer systems. Further- 
more, we provide strong evidence that in the SS phase, 
the SF condensate amplitude has a spatial modulation 
and possesses a lattice symmetry on its own. This kind 
of non-uniform condensate has been discussed related to 
solid '^He [l3| but has not been confirmed in numerical 
studies before. 

In particular, we study the 2D bilayer AF spin dimer 
XXZ model 

Hxxz = -/1^5^,, + J^Si,, .82,, (1) 

Q,2 i 

where a = 1, 2 denotes the layer index and the third term 
sums over all nearest neighbors (n.n.) in the xy plane for 
both layers. Since J >> J' (we take J' / J — 0.29 in this 
paper), the low energy physics is conveniently described 
by the interaction of the four spin states, singlet (|s)) and 
triplets S"^ = 0, ±1 (|to.±)) of each dimer. 

When A = 1, i.e., the isotropic case, the model suc- 
cessfully describes the SF phase, characterized by the in- 
plane staggered magnetization rrixy, of the spin 1/2 dimer 
compounds TlCuClg, KCuClg, and BaSiCu206 [11 El- 
The mean field (MF) ground state is 

1*0) = \{[u\s), + {-l)\ve^<^\t+),+we-^'\t_),)\{2) 



where u^ -t- u^ -I- w^ = 1 and we have neglected |io) states 
here. The order parameter is given by m^y = u{v + w) 
whose square is shown to be equal to the condensate 
density (ng, the zero momentum number density) of the 

semi-hardcore boson bl — ^~J (S'j'^j — 5^^) , which al- 
lows up to double occupancies on each site. The global 
phase 9 corresponds to the orientation of the mxy in the 
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FIG. 1: (color online) Phase diagrams of (a) the number den- 
sities Us, nt+, nto, and nt-, and (b) the order parameters no, 
S{n,n)/N, mz/2, and ml/2. The upper inset shows the de- 
tails of nt- and nto around the SS phase with enlarged scale. 
The lower inset compares the QMC (symbols) and the MF 
(solid lines) results of mJ/2 as a function of A. 



xy plane. Ignoring |t-), i.e. let w = 0, it is a hard- 
core model with n.n. repulsive interaction V and hop- 
ping t. For the case of square lattice, V ~ t — J/2. 
This Haniiltonian preserves the particle-hole symmetry, 
which is broken when |t_) is taken into account and thus 
an asymmetric phase diagram of mxy{h) is observed. 

By increasing A > 1, the repulsive interaction between 
the triplets enhances. The ratio of interaction energy and 
kinetic energy is determined by A = V/t. Therefore, a 
large value of A will lead to a solid phase where the struc- 
ture factor has a finite peak at certain wave vector. The 
compound SrCu2 (603)2 |9( has shown such solid phase 
in its magnetization plateaus. When both interaction 
energy and kinetic energy play crucial roles, the coexis- 
tence of both SF ordering and solid ordering is possible. 
Our numerical simulation indeed shows the presence of 
SS phase from A w 3 to 4. In our QMC simulation, we 
employ the stochastic series expansion approach [l5j on a 
2D bilayer square lattice for different linear sizes L = 10, 
12, 16 and 20, with temperature T inversely proportional 
to L such that at L = 20, T/J = 0.02154. 

In Fig^ we show the number densities (la) and order 
parameters (lb) as a function of magnetic field h for A = 
3.3. Let us discuss the different phases starting from low 
field h. 

superfluid I (SFI) - At low field, the phase is character- 
ized by the only non- vanishing order parameter, the stag- 
gered magnetization along z-axis m^ = {J2ii~^Ti'^ii ^ 
S2i))/N (N is the number of dimers), which is field in- 




FIG. 2: (color online) Finite size scaling of (a) S(7r, tt)/N and 
(b) no. In (c), results of S{n,n)/N (open symbols) and no 
(filled symbols) for L = 12, 16 and 20. 



dependent. This is a SF phase forming by |to) state 
in the background of \s) state with a MF wavefunction 
I'^to) = J2i{uo\s)^ + e''^(-l)Vo|io)j), where uq and vq 
are given by minimizing the MF energy and have the 



form Uq ~ i(l 
ml = 2uqVo cos (j) 



and Vq 



'u^y 



TKri "'^"^ "0 — 2^-" " iKT'- Then, 
(6q) -|- (bo), with 6q is the hardcore 
boson operator that creates a \to) state. At finite field, 
the global phase cf) is fixed to zero so that the largest 
ml aligned to the field direction. The inset of Fig^Jj 
shows an excellent agreement of MF ml/2 to the values 
obtained from QMC for different A's. Even for A = 1, 
that is relevant to many spin dimer materials, as long as 
1 > J' /J > 1/4, so Wo > 0, there is still a small but finite 
ml in this low field phase. 

Superfluid II (SFII) - When h/J > 1.87, the strong 
Zeeman energy favors |t+) instead of |io) and leads to the 
condensation of \t+). In FigQIb), the condensate density 
no = m^ of boson 6"^, jumps abruptly to about 0.2 as 
ml drops to a small value. Note that nto is small but 
still finite so that |io) behaves as impurities. This phase 
reappears at even higher fields. 

Supersolid (SS) - At h/J ^ 2.10, no begins to 
reduce rapidly while the structure factor S{Q)/N — 
^^Ani^t+nj^t+e.^^^^^) /N"^ at the wave vector (7r,7r) of 
\t+) increases. The single peak of S{n,T:)/N implies a 
checkerboard solid ordering of |i+) in the system. The 
transition from SFII to SS is also signaled by the kinks 
observed in the m^, no, rit-, and nto (FigDl- Within 
2.2 > h/J> 2.1, both SF and soHd order parameters re- 
main finite and constitute the SS phase. Finite size scal- 
ing shown in FigEl indicates that the SS phases is stable 
in the thermodynamic limit. To address the concerns of 
phase separation, i.e. the possibility of having separate 
domains of checkerboard solid and SF, we measure the 
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FIG. 3: (color online) Histograms of the averaged mz/2 of the 
whole lattice ((a) to (c)) and of a small area ((d) to (f )) of size 
4x4 (unshaded) and 6x6 (shaded). A = 3.5. For comparison, 
histograms of SFII and QS are also presented. 



average vriz within a small area of size 6x6 and 4x4 of a 
20x20 lattice. Since niz differs for different domains, the 
small area measurement should give two distinct peaks 
in the histogram in the case of phase separation. Shown 
in Fig 131 are the histograms of TOz/2 for selected fields. 
The distribution is wider in 6x6 case when compared to 
that of whole 20x20 lattice where fluctuation is averaged 
out. The single Gaussian peak at h/J = 2.15 indicates 
that there is no domains of different magnetization and 
rules out the possibility of phase separation. 

Quantum solid (QS) - Increasing the field further 
{h/J > 2.20) will half-fiU the system with 6+ so that 
rUz = 0.5. Because of strong repulsive interaction be- 
tween nearest neighbor 6+, adding or removing extra 
boson to this checkerboard solid state requires a finite 
energy. This causes the stability of the state and the 
plateau in the m^ curve. Differing from the classical 
value 5(7r, n, )/N — 1/4, the calculated structure factor 
in the solid phase is about 0.13 for A = 3.3. The m^ of 
each sublattice are found to be 0.890 and 0.110 respec- 
tively (Table P). As A increases, S{ti,tt)/N approaches 
the classical value of 1/4. At higher fields, the melting of 
the solid is of first order and the SFII phase reappears. 
nt+ keeps increasing while all other states reduces to zero 
at /i = 3.50, where all the spins are fully polarized. 

Here we address the question of non-uniform conden- 
sate in the SS phase. Table ^ shows that the m^ dif- 
fers in different sublattices A and B, and breaks the 
translation symmetry. To examine the condensate more 
closely, we define the sublattice condensates ua and ub 
as Ua = X^iGQ^^I)/-^' where a = A,B. In our simu- 
lation, we compute the products of sublattice conden- 
sate nap = X^iea jefli^I^i)/-^^- ^^ shown in Table Q] 
uaa ^ riBB and n^s = «SA = y/nAAnBB in the SS 
phase. Thus naf3 = nafip, and ua — 0.130, ub = 0.148. 
It indicates that the condensate has a checkerboard or- 
derin g on its own, which is the intrinsic nature of a SS 
state [iSj . This is different from the idea of having a SF 
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FIG. 4: (color online) Sublattice condensates uaa and ubb 
in the SS phase extrapolated to T -^ and L ^ oo. 



component "on top" of a bose solid. The checkerboard 
order of the condensate disappears in the usual superfluid 
phase SFII, while the condensate amplitude on both sub- 
lattices vanishes in the QS phase as expected. The result 
for different h is plotted in Fig. ^ The SS phase, instead 
of being a crossover between SFII and QS, is clearly a 
new phase with second order phase transitions from SFII 
and QS respectively. Nevertheless, a more detail study 
of scaling behavior is needed to determine the scaling ex- 
ponents. This result also rules out phase separation. 

Fig[5ja) and (b) show the plots of order parameters as 
a function of h for different A and the phase diagram of 
A vs. h, respectively. The isotropic case, A = 1, is simi- 
lar to those observed experimentally on other spin dimer 
compounds. Enhancing the anisotropy to A = 3, there 
is a small drop in no around m^ = 0.5 while S{'k, '!t)/N 
remains zero. For even larger repulsive interaction at 
A — 3.1, the drop develops further with a finite peak of 
S{TT,Tr)/N at rriz = 0.5. Note that it is a commensurate 
SS that is made possible by the semi-hard core nature 
of the spin bosons. The finite peak of S{TT,Tr)/N finally 
develops to a plateau together with the ruz in the QS 
phase. As a consequence of lowering energy of \to) state 
with increasing A, the phase boundary of SFI to SFII 
moves to the higher fields and shrinks the region of SFII 
and SS which eventually vanishes at A ~ 4 and leaves a 
first order transition from SFI to QS phase. As shown 
in FiglSfb), the SS phase is stable in a parameter region 
with 4 > A > 3. 

The next question is whether these phases can be re- 



SFII 
SS 
QS 



4/2 



0.224(5) 
0.397(1) 
0.4449(2) 



nif/2 



0.225(5) 
0.0831(9) 
0.0549(2) 



riAA 



0.0384(7) 
0.0169(5) 
0.0017(1) 



riBB 



0.0384(6) 
0.0219(5) 
0.0036(1) 



nAB 



0.0384(7) 
0.0195(5) 
0.0026(1) 



TiBA 



0.0387(6) 
0.0194(4) 
0.0027(1) 



TABLE I: Sublattice magnetizations and condensates for SFII 
{h/J=2m), SS [h/J =2.14) and QS {h/J = 2.24) phases. 
A = 3.3 and L=20. The small values of n^/s in QS phase is 
due to the finite size effect. 
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FIG. 5; (color online) (a) Order parameters as a function of 
h/J for different A with lattice size 20x20x2. The inset shows 
the isotropic case (A = 1). (b) Ground state phase diagram 
(extrapolated to i — > cxa) of A vs. h/ J . Solid (dotted) lines 
denote the first (second) order phase boundaries. 



is approximately given by 



Si,iS2j- ~ Si^Sij — t:{ShSi , + SuSi 



(3) 



The inclusion of J" now breaks the spin rotational sym- 
metry and leads to an XXZ model with effective n.n. 
coupling of J* = J' — J" and anisotropy A given by 
J* A = J' + J". Therefore, if J" = J' /2 one has a spin 
dimer XXZ model with A = 3, close to the SS phase. 

In summary, we demonstrate that the spin dimer XXZ 
model, a natural semi-hardcore boson system with de- 
fects, contains a SS phase that characterized by the bi- 
parte condensate density. The anisotropy can be a con- 
sequence of including the n.n.n. coupling among dimers. 
While it is suggested 7] that SS state of hardcore bosons 
can be tested in triangular optical lattices, we propose 
that spin dimer compounds may be a natural place to 
realize the SS state of semi-hardcore bosons. 

We acknowledges financial support by the NSC 
(R.O.C.), grant no. NSC 94-2112-M-001-003. 



alized in real material. There are two issues, however, 
that have to be considered. The first is the anisotropy 
that introduced in Hxxz only applies to the interdimer 
coupling but not the intradimer coupling. This kind of 
selective anisotropy is hard to realize. The second issue 
concerns the large value of A ~ 3 for the SS to be sta- 
ble. In usual compounds the XXZ anisotropy is caused 
by the spin-orbit coupling, which couples the spin to the 
crystal structure, and the effective A is much smaller 
than 3. To the first issue, one may instead consider a 
S — I compound with a single ion anisotropy D^{S^)^, 
where D plays the role of intra-dimer coupling J and 
the states S^ = ±1 and S^ = plays the role of |i±), 
and |s), respectively. There is no corresponding states 
of \to)- Therefore, SFI is replaced by 5*^ = state on 
each site. Together with spin-orbit coupling that real- 
izes the XXZ anisotropy, one has essentially the same 
phase diagram as in the case of Hxxz, although a large 
enough anisotropy is still needed. On the other hand, 
there is another route to look for the realization of the SS 
phase. Consider two neighboring i_|_ states, the inclusion 
of n.n.n. coupling J"Si_iS2j {i,j are n.n.) will simply 
enhance the repulsion between \t+) states. Contrastingly, 
if |to) is ignored, J" reduces the hopping between the 
neighboring |i_|_) and \s) states due to the antisymmetry 
of \s). In bond operator representation, when the |io) 
states is ignored, spin operators 5*^ 2 ~ k(^+^+ ^ i-t-) 



[9 

[lo; 
[11 

[12: 
[is; 
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[is; 



and S'+2 « ±^(stt 



V2^ 



tls) 



Then the n.n.n. coupling 
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